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Abstract 

Based on symmetries Taub-NUT shares with Bertrand spacetime, we cast it as the latter 
with magnetic fields. Its nature as a Bianchi-IX gravitational instanton and other related 
geometrical properties are reviewed. We provide an easy derivation and comparison between 
the spatial Killing-Yano tensors deduced from first-integrals and the corresponding hyperkahler 
structures and finally verify the existence of a graded Lie-algebra structure via Schouten- 
Nijenhuis brackets. 

Contents 

1 Introduction 2 

2 Conserved Quantities 3 

2.1 A dynamical-systems description of Taub-NUT. 5 

2.2 Holten Algorithm description. 7 

2.2.1 Some basic Killing Tensors . 7 

2.2.2 Holten algorithm for Taub-NUT . 8 

3 Bertrand spacetime dualities 9 

3.1 Bertrand spacetimes with magnetic fields . 9 

3.2 Kepler-Oscillator duality. 11 

4 A review of geometric properties 12 

4.1 Taub-NUT as a Darboux-Halphen system. 13 

4.2 Curvature and anti-self duality. 14 

4.3 Topological Invariants. 16 

5 Killing-Yano tensors and the Taub-NUT metric 17 

5.1 Yano and Stackel tensors . 17 

5.2 Euclidean Taub-NUT . 18 

5.3 Graded Lie-algebra via Schouten-Nijenhuis Brackets. 19 

6 Hyperkahler structure and the KY tensors 20 

7 Discussion 22 

8 APPENDIX 24 

8.1 Basic Killing tensors from Holten’s Algorithm. 24 

8.2 The Bohlin transformation. 25 

8.3 Double derivative of Killing-Yano tensors . 26 

















1 Introduction 


The Taub-NUT [1] is an exact solution of Einstein’s equations, found by Abraham Huskel Taub 
(1951), and extended to a larger manifold by E. Newman, T. Unti and L. Tamburino (1963). It 
is a gravitational anti-instanton with corresponding SU (2) gauge fields, frequently studied for 
its geodesics which approximately describe the motion of well seperated monopole-monopole 
interactions. As a dynamical system it exhibits spherically symmetry, with geodesics admitting 
Kepler-type symmetry, implying first-integrals such as the angular momentum and Runge-Lenz 
vectors respectively. Witten’s prescription [2] realized Taub-NUT space as a hyper-Kahler 
quotient using T-duality. This construction has a natural interpretation in terms of D-branes 
[3], serving as an important example in string theory. 

The Bertrand spacetime metric, formulated by V. Perlick [4] is also spherically symmetric 

ds 2 = h(p) 2 dp 2 + p 2 (dO 2 + sin 2 9 dcf) 2 ) — ^ (1.1) 

derived from Bertrand’s Theorem, describing stable and closed geodesics with periodic 
orbits. Upon comparison, Euclidean Bertrand spaces and Taub-NUT spaces, appear quite 
similar apart from magentic monopole and dipole interaction of the Taub-NUT. This implies 
dynamical similarities, manifested through similar first-integrals characterizing their motion. 
It also implies that Taub-NUT possibly exhibits Kepler-Hooke configuration duality. 

Consequently, we try to find first-integrals similar to those associated with central-force 
motion under potentials involved in Bertrand’s Theorem: the angular momentum and Laplace- 
Runge-Lenz vector. Since we are interested in the dynamical aspects of Taub-NUT spaces, 
our attention is directed toward geodesics and Killing tensors. Naturally, we will be looking 
at Killing tensors affiliated with Runge-Lenz-like vector. They obey the equation: 

V (a ^i) 62 ...6„ = 0 (!- 2 ) 

Such tensors are the Killing-Stackel tensors which are symmetric under index permutation 
and the Killing-Yano tensor. The Killing-Yano tensors are antisymmetric under index permu¬ 
tation, and their square gives the Stackel tensor, like the antisymmetric tensor whose square 
gives the Runge-Lenz-like quantity as we shall see. Such Killing tensors exhibit quaternionic 
algebra, implying a connection to Hyperkahler structures associated with the metric. 


We start in section 2, with preliminaries on mechanical systems with magnetic field in¬ 
teractions, then compute first-integrals similar to the angular momentum and the Laplace- 
Runge-Lenz vector, in forms specifically for the Taub-NUT. We deduce such first-integrals 
using equations of motion and analytically using a momentum polynomial expansion. 

In section 3, we compare Taub-NUT metric to Euclidean Bertrand spacetime with magnetic 
monopoles and dipoles. Demonstrating such a similarity allows the intensely studied Bertrand 
spacetimes to share many important properties, and conversely extend properties of the Taub- 
NUT to Bertrand spaces with magnetic fields. This helps us identify symmetries and conserved 
quantities of Taub NUT and employ its curvature properties for Bertrand spacetimes. The 
last subsection covers the conserved quantity called the Fradkin tensor under Bohlin-Arnold- 
Vassiliev transformation which are bound to have such Killing tensors embedded. 

In section 4 we derive the Taub-NUT from a special case of self-dual Bianchi-IX met¬ 
ric described by the classical Darboux-Halphen system. Then we geometrically analyze it, 
computing curvature and confirming its self-duality as a gravitational instanton. 


Rpupa = R 


X'-ypa 


(1.3) 
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This analysis helps us explore the metric as an integrable system. Finally, we will compute 
topological invariants shared with comparable Bertrand spacetimes with magnetic fields. 

In section 5, after a short introduction to Killing Stackel tensor and Yano tensors, we 
will focus on the latter. After a brief overview of their properties, we will attempt to find 
them embedded within conserved quantities. Then, we see if it exhibits a graded Lie-algebra 
structure that decides if higher order Killing-Yano tensors can be constructed from it. 

Finally, in section 6 we derive hyperkahler structures of the Taub-NUT. Then we compare 
them to the Killing-Yano tensors to see if they also exhibit quaternionic algebra. 

In the last section we conclude our works and discuss the possibilities of further research 
along the line of the present article. The appendix contains detailed computation regarding 
Killing tensors following Holtens algorithm, a brief review of Bohlin’s transformation of the 
Taub-NUT and a double derivative expansion of Killing Yano tensors. 


2 Conserved Quantities 

In classical mechanics it is important to identify constants of motion called conserved quantities 
or first-integrals of the system. In the theory of integrable systems, all first-integrals are in 
involution or commute with each other within the Poisson Brackets, with at least one integral 
definitely being available. Given a n + 1-spacetime metric with t as cyclic variable: 

ds 2 = gij(x)dx l dx^ + gu{x)dt 2 + 2 gi t (x)dx l dt (2.1) 

parameterized as t = r, we will have the Lagrangian and a conserved quantity q: 

h =i = + ^g tt (x) + g it {x)x l q = 00 = 9tt + 9ux l (2.2) 

The Hamiltonian is given by the Legendre transform H = Yhk^t ’M%x k ~ A so that: 

H = i g lJ (x)(pi - g it (x)) (pj - g jt {x)) - \gtt(x) p i = (2.3) 

In Hamiltonian dynamics, a conserved quantity Q commutes with the Hamiltonian H, a 
first integral resulting from timetranslation invariance, within the Poisson Brackets: 

{Q,H}= 0 (2.4) 

However, this prescription is not gauge covariant for systems with gauge interactions. To 
better understand why, consider the following metric with scalar potential U(x): 

ds 2 = 5ijdx l dxi — — ~-^-^—^-dt 2 (2.5) 

m 

where t is cyclical. Under the parameterization t = r, the lagrangian, Hamiltonian and 
Hamilton’s dynamical equations for particles in presence of scalar potentials is given by: 


m 


L i =l = y* 2 ~U(x) 


1 


H = ^—p 2 + U(x) => < 

2m 


. _ dH _ p 

Op m 
OH 

P = -- = -W(*) 


( 2 . 6 ) 

(2.7) 
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For this system without magnetic fields, the fundamental brackets are: 

{x l ,pj} = 5 lJ {x\x 3 } = {pi,pj } = 0 (2.8) 

Now, for charged particles in U(l) gauge fields from magnetic dipoles alone, without scalar 
potential, the metric is: 


ds 2 = 8ijdx l dx 3 — (df 2 — 2Ak(x)dx k dt) 


m 


(2.9) 


so the corresponding Lagrangian and Hamiltonian for t = 1 are given by: 


Li = i = A ( mx 2 — i 2 + 2 A(x).xi 


q=\ 


fdL 


\dtj i= i 


= 1 — A(x).x 


H 


2m 


( p-A{x)Y 


( 2 . 10 ) 


For charged particles in the presence of magnetic monopole and dipole U(l) gauge fields 
without scalar potential, the metric is: 


ds 2 = 5ijdx l dx 3 —— (dt — AAx)dx k ) 2 
m ' 

so the corresponding Lagrangian and Hamiltonian for t = 1 are given by: 


( 2 . 11 ) 


h =i = 


mx 


— (l — A(x).xy 


q=\ 


fdL 


\di )i =l 


= 1 — A(x).x 


( 2 . 12 ) 


h = ^-qMx))- 


Now let us consider a Kaluza-Klein modihcation of this spacetime, such that we include 
another cyclical co-ordinate if that is periodic along with magnetic field components coupled 
with it. This would result in a 4 + 1 spacetime from a 3 + 1 one given by: 


ds 2 = 5ijdx l dx 3 H- (dif + Ak(x)dx k ) — (l + 2 V(x))dt 2 (2-13) 

so the Lagrangian and Hamiltonian for t = 1, ignoring constant additive terms are: 

1 r i c) L 

mx 2 + (ip + A(x).x ) —V(r) q = —- = ip — A(x) 

df' (2.14) 


L = 


2 L 


H = ^(p - qA(x)) 2 + V(r) 

where q is a conserved charge. The corresponding Hamilton’s equations are: 
d H p — qA 


x = 


dp 


m 


dH q , w . 

p = ~— = -(VA).(p-qA) -VV 


dx 


m 


(2.15) 


Since the potentials are gauge dependent (A —> A + VA), the momenta therefore must be 
so as well (p —> p + q'V A). Then, we must write gauge invariant momenta and express the 
Hamiltonian in its gauge invariant form. 


n 2 

H=—+V(r) 


II = p — qA 


(2.16) 
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Any function and partial derivative operators in gauge invariant forms can be written as: 


/(*.p) - 

-»• f( x 

,n) 


d 

dW 

d 

d 

L 

dx 1 

dx 1 

dW 

dx 1 

d 

dW 

d 

d 

dp 1 

dp 1 

dW 

dp 1 


= -qdiA 

d 

~ dw 


d 


+ 


d 


2 dW ' dx 1 


(2.17) 


No explicit dependence on p 


with which the fundamental brackets become: 


11 ; } 6 ij {x\x j }=0, 


(2.18) 


{llj, Ilj} = -qFij 

where it is interesting to note that the new Poisson Brackets between the gauge covariant 
momenta are non-zero, as opposed to the usual case. This is a classical analogue of Ricci- 
identity (in the absence of torsion). We can furthermore redefine the Poisson Brackets as: 


r, 1 = ^L . .^L _ p d f dg 

1 a™ ‘ ott ;itt ’ nrr " ;itt 


(2.19) 


dx on <9n dx 1 J on sn 

Now that we have redefined the Poisson Brackets to make Hamiltonian dynamics manifestly 
gauge invariant in the modified bracket, we can proceed to analyze the conserved quantities 
in a general gauge invariant form. This is done by the Holten Algorithm as shown in [5] and 
[6] discussed later as we shall see. 


2.1 A dynamical-systems description of Taub-NUT 

The Euclidean Taub-NUT metric as shown in [1] is given by: 

ds 2 = f(r){dr 2 + r 2 (d0 2 + sin 2 6 d(j ) 2 )} + g(r)(di/; + cos 9d(fi)~ 

4 M (4 M) 2 

where f(r) = 1 + — g(r) = 


( 2 . 1 . 1 ) 


ds 2 

For later reference, taking ds 2 = - we shall re-write the above metric into this form : 

’ B 4 M 


ds 2 = V(r) 5ij dx l dxi + V 1 (r)(d-0 + A.dx ) 

where V(r) = —— -|— A.dx = cos Odd) 

v ’ 4 M r 


( 2 . 1 . 2 ) 


We now consider the geodesic flows of the generalized Taub-NUT metric given by (2.1.1), 
for which we can compose the Lagrangian: 


C = i f(r){r 2 + r 2 [0 2 + sin 2 0 <j > 2 )} + ^g(r)(ip + cos 0^»)“ (2.1.3) 

We can further re-write the Lagrangian (2.1.3) into 3-dimensional form with a potential, 
as in (2.14), independent of the ijj as: 

£ = ^f(r)\x\ 2 + ^g(r)(ip + A.x) 2 - U(r) (2.1.4) 

where the momentum can be written as: 


P = 


dC 

dx 


f(r)x + qA 


II = f(r)x = p — qA 


(2.1.5) 
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Spaces with the metric ( 2 . 1 . 1 ) exhibit SU( 2 ) x U(l) isometry group. Given that we have 
at least 2 cyclical variables ip and (p, we will have the following 4 Killing vectors given by: 

D 0 = dq, 


COS <p 


D i = — sin <p do — cos (p cot 6 dq, H- 

sm 0 

sin^ 

D 2 = cos 4> dg — sin q i cot 6 cb, H—:—- 

sm (J 


d , 




( 2 . 1 . 6 ) 


<9,/ 


^ 3 = df 

where Ho commutes with all other killing vectors, while D\, D 2 , exhibit the SU(2) Lie 
algebra given by = — £ij k D *.. Since ip is cyclic, we have a conserved quantity: 


d£ 

q = —- = g(r) (ip + cos 0 (p) = g(r) (ip + A. sc) = const 
dip 


(2.1.7) 


known as the relative electric charge. The symplectic 2-form u> and energy 6 for the Taub- 
NUT system in Cartesian co-ordinates are: 

3 3 


uj = - (ujq + qF(x)) , k dx^ A dx k = E d(pi — qAi(x )) A dx l = E cTR A dx l 

i=l i=1 

3 

= ^ dpi A dx l — EijkX 1 dad A dx k (2.1.8) 

i,j,k 


i= 1 


w = + Aw W = - E ^ (2-1.9) 

k 


2 /(r) 2 g(r) 


Consequently, the Hamilton’s equations are given by: 

fc = ^ ri = {n,ft} fl = a(r)* + 


2 : 


r 3j( r ) 


* x n — W(r) 


where 


a(r) = 




|n | 2 + 


0(0 


( 2 . 1 . 10 ) 


2 (/( 0 ) 2 ( ff (r )) 2 

Using these equations, we find angular momentum in presence of magnetic fields to be: 


dx „ dll \ q 

— x II + sc x —— = — 

dt dt J r 3 f(r ) 


[sc x (sc x II)] = g 


(sc.sc) sc x 


d (x 


d ( „ x , 

— (x x n + 5 - 1 = 0 


J = x x II + g— 


= 


x 


dt\r 


( 2 . 1 . 11 ) 


The cyclic variable allows reduction of the geodesic flow on T(M 4 — {0}) to a system on 
T(R 3 — {0}). The reduced system’s rotational invariance implies it must have a conserved 
energy, angular momentum and vector K analogous to the Laplace-Runge-Lenz vector: 

- 2 x 1 n 2 


1 n 2 

n = --— + 


2 f(r) \2 g(r) 


1 q* 


+ U(r) = 


2/(0 


+ W(r) 


x 


J = x x II + q— 


1 


K = -K^x" = n x J + ( T 4 mE ) - 


Q 


4 m 


x 


( 2 . 1 . 12 ) 

(2.1.13) 

(2.1.14) 


This concludes the detailing of conserved quantities of the Taub-NUT from a dynamical 
systems perspective. Now we shall proceed to consider a systematic analytic process that 
describes conserved quantities as power series expansions of momenta. 
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2.2 Holten Algorithm description 


One way of analytically obtaining conserved quantities that are polynomials in momenta is by 
writing them in a power series expansion involving the gauge invariant momenta: 

Q = + C'f ) (r)tf + icf (r)iriF + ^C$(r)rf WU k + .... (2.2.1) 

where all the coefficients of momenta power series are symmetric under index permutation. 
Applying this to eq (2.4), we can obtain the relations for each coefficient by matching the 
appropriate product series of momenta for both the terms. 


II 

5 

EK r 

f[n‘,nJ}n, + {c<”>n n ' 

% V(r)} 

= 


n 

w « 



n n ' - lOwf R’ + KjV(r) 

n 

n 


{m} 


({ ™},k+ 

0 

The equations we will 

(i) 

get up to the 3rd order setting = 

0 V i>3 

order 0: 

0 = C$d m {V{r)) 



order 1: 

ViC® 

= qF ij cf ) +c[fd j {y{r)) 



order 2: 

v.cj 11 

+ VjCV = q{F im C^ m ] + 

F C ^ 

1 jm^ 

m r) 

order 3: 


+ V*C f + VjClf = 0 



Now we will turn our 

attention to 

some familiar conserved quantities. 



( 2 . 2 . 2 ) 


(2.2.3) 


2.2.1 Some basic Killing Tensors 

Using the above relations for the various terms, we will now look at some familiar ones that 
we have already studied in classical mechanics. 


Angular Momentum 

The conserved quantity that results from the 1st order term of the Holten series alone is: 

qW = cf } ir = -g im (x)e m ok e k xm l 
=> L.e = -(e o - fc nv)0 fc = (x x n ).0 

L = x x n (2.2.4) 

This eventually becomes the conserved quantity known as the angular momentum. 


Laplace-Runge-Lenz vector 

On the other hand, the conserved quantity from the 2nd order term of the series alone is: 

Q (2) = = {|n| 2 (n.®) - (n.*)(n.n)} 

=> N.n = {|ll| 2 a: — (n.;c)n}.n. = {n x [x x Il)}.n 

AT" = II x (x x n) (2.2.5) 

This quantity is a term contained in another conserved quantity known as the Laplace- 
Runge-Lenz vector. Having found the two familiar types of conserved quantities, we can now 
proceed to see what it looks like for the Taub-NUT metric. 
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2.2.2 Holten algorithm for Taub-NUT 


Now, for the Taub-NUT metric, we have (2.1.12) giving the Hamiltonian. This can be written 
in dimensionally reduced form as: 


U = - i \U\ 2 + f{r)W{r) 


w ^ = u ^ + W) + W)- £ 


( 2 . 2 . 6 ) 


(2) (3) 

From this Hamiltonian, after setting all higher orders C\- = C\,- k = 0, we get the modified 


1 st and 2nd order equations to be the following: 

order 1: d;C (0) = qF t] C (l)j 

order 2: V i cf ) + VjCf } = 0 

The constraint equation of the 2nd order of (2.2.7) gives us: 

Cf } = g im (x)e m jk 0ix k 

diC (Q) = \eijk e j nm x k d m x n = ^[ix(9xj)]. = -| [r 2 0 - (x.O) x \, 


ViCW =q{ 


r 1 

(x.e)xi 


C® = qdi~— 


ry* ryuJ J ry* 

Thus, we have the overall solution, and the corresponding conserved quantity: 

Q = J k e k = c ( °) + cf } it = ( - g im (x)s m jk x j U l + q^j d k 


X 


J.6 = [ x x II + q— .6 


J = x x II + q- 


x 


(2.2.7) 

( 2 . 2 . 8 ) 

(2.2.9) 

( 2 . 2 . 10 ) 

( 2 . 2 . 11 ) 


( 2 ) 

However, if we explore upto the 2nd order, setting C - ^ 0, we will return to the equations 

( 2 ) 

in (2.2.3). For the 3rd order equation, the solution for C\- is given by (8.1.7), so that: 


c\f — (2 gij[x)nk - gik(x)rtj - g k j(x)ni)x k 
= g(F it cff + 


Eventually the other co-efficients are given by: 

f(l) _ TF /"((2) 


( 2 . 2 . 12 ) 

(2.2.13) 


P P<( 2 ) _ _ 0 _ . („ rr m \ _1_ C ■ _1_ C ggk n ^3 

-Tik^kj Atijn 3 y lj m Jy ) 1 c-ikn 3 ,l j 1 c- ikn 11 ^ 3 


rj* k 


X i 

3 


( nxx)i 


VjC-^ + VjC^ = q\ ~^(n x x)i + ^(n x aj)j 


,(i) 


x 


Here, one can choose to insert extra terms: 


v^ + v.cr = -</<! V. 


( 1 ) _ 


'3 ' J 3^i - j I J l 


/ Ei km n k x rn 


+ Vj 


X C~ ■ 1 ry-fk rpTTl 

I c jkm ,L ^ 


Thus, we can easily see which term on the RHS corresponds to what on the LHS, allowing 
us to solve for the 1st order and zeroth order coefficients from (2.2.3) : 

Cf } = -lg im (x)e m jk n k x i 


(2.2.14) 









v,c<»> = -„ 2 rt - 


(; X.njXi 


+ {2(n.x)8i j -niXj-x i n j }dj ( f (r)U(r) + q 2 £-Sf(r) 


In the case of the generalised Taub-NUT metric, the most general potentials admitting a 
Runge-Lenz vector are of the form: 


U(r) = 


1 


q 


f(r ) \2r 2 r 


+ - + !)- 


q 


2. g(r ; 


+ £ 


ViC^ =j8 -- 

r 


( n.x)xi 


C<°> = Pm— 

r 


For integrability, we require the commutation relation: 


[di,dj\cW =0 


A ( f(r)W(r) — 


q 2 g 2 

2 r 2 


= 0 


ft It \ P , 

f(r)W(r)--^- T = -+ 7 


ft MT 71 \ q2g2 , P 1 
f(r)W(r) = -^r + ~+l 


Thus, this overall conserved quantity is given as: 

Q = R k 0 k = c (0) + cf ] ir i + d 2) u l IF 


R.n = (n x (a; x n) — -x x U + /3 


P,1 G 


.n 


R = n x J + (3- 




(2.2.15) 

(2.2.16) 


(2.2.17) 


(2.2.18) 


Now we will take a detour to look at some details regarding the Runge-Lenz vector. 


3 Bertrand spacetime dualities 

In Newtonian mechanics, there are only two potentials allowing stable, closed and periodic 
orbits: Hooke’s Oscillator (V(r) = ar 2 + b ), and Kepler’s Orbital Motion (T(r) = - + b) 
potentials. There is a relativistic analogue, given by the corresponding metrics in [4], describing 
spherically symmetric and static spacetime, with bounded and periodic trajectories. The Taub- 
NUT is one example of a spherically symmetric spacetime. Naturally, one would ask how it 
compares with the Euclidean Bertrand spacetime (BST) metric with magnetic fields. 


3.1 Bertrand spacetimes with magnetic fields 

The Bertrand spacetime metric is given by (1.1). If we take the Euclidean version and include 
magnetic monopole and dipole interaction terms, then the metric becomes like (2.11) as: 


ds 2 = h(p) 2 dp 2 + p 2 (d0 2 + sin 2 6 dtp 2 ) + (dt + Aidx 1 ) 2 

T{p) 


(3.1.1) 


If we recall, the Taub-NUT metric was given by (2.1.1). To see how they are comparable, 
we shall attempt a co-ordinate map. 

1 


fir) dr = h{p) 2 dp- f(r)r 2 = p 2 g(r) = 


r (p) 


t = if} + k 


dr h{p) dp 
r io 


f dp AH 

r = r o e J p 


9 













Thus, we can suppose that Taub-NUT metric resembles Bertrand spacetime with magnetic 
fields. We can also proceed the other way around, starting with the generalized Taub-NUT 
metric and proceeding toward Bertrand spacetimes by applying appropriate potential power 
laws shown in [7]. Thus, like the BSTs, there are two Taub-NUT configurations: 


1. Hooke’s Oscillator configuration: 


2 . 


fo(r) = 

Kepler’s orbital configuration: 

f/c{r) = 



9o(r) 

r 2 (ar 2 + b) 
cr 4 + dr 2 + 1 

ar 2 + b 

a + br 

9K.(r) 

r (a + br) 

r 

cr 2 + dr + 1 


Evidence for the duality between these two configurations of the metric can be clearly 
demonstrated. To study Taub-NUT space duality, we confine motion to a cone (6 = const). 
This is permissible because of the conserved angular momentum (2.1.11), for which [1, 7, 8] 


J.e r = | j| cos 0 = const => 6 = const (3.1.2) 

This allows us to reduce the problem to 2-dimensions by rendering 6 a constant co-ordinate, 
allowing us to write the metric as: 

ds 2 = f(r)(dr 2 + r 2 a 2 dcp 2 ') + g(r)(dip + pi d<p) 2 a = sin 9, (3 = cos 6 (3.1.3) 


We shall represent the co-ordinates as Z = x + iy, £ = X + iY, where \Z\ = rcos| and 
perform Bohlin’s transformation [9] of the Oscillator metric (Z —> £ = Z 2 ). This requires the 
complex co-ordinates defined for self-dual Euclidean spaces [10], where (9 = const ) (3.1.2): 


Z = x + iy = \Z\ exp 


; (ip + 4>) 


Z —>• £ = Z 2 = |Z| 2 exp \i(il> + (/>)] 


£ 


X + %Y = |£| exp 


;(* + *) 


<h = 2 (j), -0 —>- x = 2g/> 


(3.1.4) 

(3.1.5) 


(ds 2 ) 0 


(a\Z\ 2 + b)\dZ\ 2 + 


\Z\ 2 {a\Z\ 2 + b) 
c\Z\ A + d\Z\ 2 + l 


(dip + /3 dtp ) 2 


(3.1.6) 


W +<o m 2 ++ f> #)’ 


Z ->£ = z 2 

m 2 + 


1 f a l£l + & l^|2 


Ifl 


4> —>■ T = 2 <p, ip —»■ x = 2-0 

ICIHCI + &) /-j , o 
cK |2 +rf | e | + 1 (^ + <3 *) 


Then we can compare with the Kepler system in presence of magnetic fields: 


(ds 2 ) K 


b\Z\ +a 
\Z\ 


\dZ\ 2 + 


\Z\ (b\Z\ + a) 
c\Z\ 2 + d\Z\ + 1 


(dx + pd^f 


(3.1.7) 


showing that aside from a factor of —, a variable swap a b completes the transformation, 


and thus, the two configurations of Taub-NUT are also related via Bohlin’s transformation like 
Bertrand spacetime. For various settings of the constants, one can get different configurations 
of spacetime, as shall be described in the following table. 


10 
















Systems for various settings 


Type 

a 

b 

c 

d 

f(r) 

g{r) 

System Name 

1C 

0 

1 

1 

9 

1 

r 2 

MIC-Zwangier 


(1 - r) 2 

1C 

0 

1 

0 

2k 

1 

r 2 

MIC-Kepler 

q 2 

1 - 

y 

o 

0 

1 

k 

0 

1 

r 2 

MIC-Oscillator 

q 2 

1 + 4r 4 

qZ 

1C 

4m 

1 

0 

1 

4m + r 

(4 m) 2 r 

Euclidean Taub-NUT 

4m 

r 

4m + r 


1C - Kepler, O - Oscillator 


3.2 Kepler-Oscillator duality 

In the study of central force problem, we learn that the Kepler and Oscillator systems are dual 
to each other according a duality map demonstrated in [11] and [12]. This is summed up in 
Bertrand’s theorem which describes them as the only systems with stable, closed and periodic 
orbits. Thus, curved Bertrand space-times are classified as Type I and Type II, representing 
Kepler and Oscillator systems respectively. 

If we start with the 2-dimensional simple harmonic oscillator described by x l = —w 2 x®, we 
are reminded of a conserved tensorial quantity, known as the Fradkin tensor: 

T ij = pV + kxV i,j = 1,2 (3.2.1) 

Any conserved quantity can be obtained by contracting the Fradkin tensor over its two 
indices by any chosen structure, ie. 

Q = MijT ij (3.2.2) 

This quantity is symmetric under index permutation. Its complex counterpart is given by: 

r ^'z a z b = G^zbTij z = (3.2.3) 



According to the Arnold-Vasiliev duality [13], a co-ordinate transformation and re-parametrization 
of the first two complex Fradkin tensors will give us the Laplace-Runge-Lenz vector. 

A = pxL + /3- (3.2.5) 

r 
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In tensorial form, this is written as follows: 

Ai = £i ki £ l jmpkp m + = x J | (SijS km - S ik 5 jm )p k p m + ^5ij | (3.2.6) 

showing that the 1st term can be expressed in a form quadratic in mommenta. Since it 
is essentially a linear combination of Fradkin tensor components, we would prefer it to be 
symmetric in the momentum indices like its oscillator counterpart. Thus, we can write 

Ai = x J | ^ ( 25ijd km - S ik S jm - d im d jk )p k p m + j (3.2.7) 

Hence, to describe this conserved quantity of the Kepler system, we need tensors that are: 

1. quadratic in momenta 

2. symmetric under index permutation 

3. conserved along geodesics 

Our next step will be to explore such tensors in the next section. 


4 A review of geometric properties 

An instanton or pseudo-particle is a concept in mathematical physics that describes solutions 
to equations of motion of classical field theory on a Euclidean spacetime. The first such 
solutions discovered were found to be localized in spacetime, hence, the name instanton or 
pseudoparticle. Instantons are important in quantum field theory because: 

1. They are leading quantum corrections to classical motion equations in the path integral 

2. They are useful for studying tunneling behaviour in systems like the Yang-Mills theory 

Since we are considering the Taub-NUT metric defined on a 4-dinrensional Euclidean plane, 
it is worthwhile to verify if it is an instanton as well. In this section, we will proceed to analyze 
its geometrical properties exhaustively, verify if the Taub-NUT is an instanton from the cur¬ 
vature components computed from the metric, and also take a look at its topological properties. 

A variable transformation m = 2 M and p = r — 2 M of (2.1.1) gives the Taub-NUT as: 

ds 2 = — m dr 2 + 4m 2 -—— (di/j + cos 9 dcj)) 2 + (r 2 — m 2 ) ( dO 2 + sin 2 9 d<p 2 ) (4.1) 

r — m r + m y 

This can be further recast into the form: 

ds 2 = r -^dr 2 + 4 m 2r -^a 2 + (r 2 - m 2 ) (a 2 + a 2 ) (4.2) 

r-m r + m v ;v 

where the variables <7j are essentially solid angle elements in 4-dimensional Euclidean space 
obeying the following structure equation: 

da 1 = — £ l j k cr 7 A a k a 1 = —-^rf filJ x ll dx l/ (4.3) 

Now that we have identified the vierbeins, we will proceed to implement Cartan’s method 
of computing spin connections and the Riemann curvature components. Embedded within 
them are the SU(2) gauge fields and their corresponding field strengths as we shall see. 
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4.1 Taub-NUT as a Darboux-Halphen system 

The Taub-NUT is a special case of self-dual Bianchi-IX metrics [14], which are characterized 
by the classical Darboux-Halphen system. The self-dual metric and its characteristic system 
of equations are given by: 

ds 2 = (tlitl 2 tl 3 )dr 2 + ^ 2 ^ 3 (o~i ) 2 + ^-^(cr 2 ) 2 + ^^(o - 3 ) 2 (4.1.1) 

'*3 

til = 0 2 ^3 — tli (tl 2 + tl 3 ) 

ti' 2 = ti 3 rii - n 2 (n 3 + tii) ()' = ^() (4.1.2) 

tig = DlD 2 — tl 3 (tl 3 + tl 2 ) 

where tlj are parameters defined to re-write the Bianchi-IX metric into the above form 
(4.1.1) to write self-dual equations. One particular first integral of the this system [15] is: 


Q = 


(tli ) 2 


+ 


(n 2 ) s 


+ 


Os ) 5 


(ti 3 — tii) (tii - ti 2 ) (tii — ti 2 )(ti 2 — ti 3 ) (ti 2 — n 3 )(n 3 — tii) 


(4.1.3) 


In case of the Taub-NUT, we need to set tl 2 = tl 3 = tl 7 ^ tii = A. This way, we will get 
the following metric, system of equations and first integral: 


O 2 


ds 2 = tl 2 A dr 2 + A[(<t 2 )' + (cr 3 ) ] + -^-( 01 ) 


(4.1.4) 


dA 

dr 


- 2A) 



(4.1.5) 


lim Q 
n 2 — 


=a 


A2 + D_[ , im (JAZ _ShlTl] 

(A-ti) 2 A - ti |_!U^3=^ vtl 2 - tl 3 ti 2 -ti 3 y 
A 2 2tl _ / tl \ 2 

(A - tl) 2 + A — tl “ " VA-tl ) 


(4.1.6) 


Rescaling the radius and solving (4.1.5) with suitable constants of integration gives us: 


dr 

<m 1 

d (A\_ 

1 dtl 


2 mtl 2 

dr 2m 

dr \Q 2 ) 

n 2 dr 



r — m 

2 2 
r — m 



tl = 

2m 

Am 2 


(4.1.7) 


and rescaling the metric as ds = —— we get the Taub-NUT (4.2) and conserved quantity: 

ds 2 = 1 +m dr 2 + 4m 2 -—— (cti) 2 + ( r 2 - m 2 ) [(cr 2 ) 2 + (cr 3 ) 2 ] 
r — m r + m L ' ' x J 


lim Q 

O2 —1CZ3— 


= -1 - 


- fii=A 


tl 


A-tl 


r 2 — 2 mr + 5 m? 
(■r — m) 2 


(4.1.8) 


This concludes another possible symmetry of the Taub-NUT as a member of Bianchi-IX 
metrics or solutions to Darboux-Halphen systems. 
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4.2 Curvature and anti-self duality 

Now that we have identified the individual vierbeins, we shall proceed to compute the spin 
connections. We can describe the vierbeins as: 


e° = co(r)dr e l = Ci(r)a l , z = 1,2,3 (4.2.1) 

Obviously, e° produces no connection terms (de° = 0). Under torsion-free condition the 
1st Cartan structure equation (de l = —co l j A e J ) gives us the spin connections. 


i dr Cj j 

CJ o = - c 

Co 


C0 


= — £ 


jk~ 


+ I J2 2 

C i + Cj c k 

2 C} Cj 


a 


(4.2.2) 


The elaborate form of the spin connections is used to keep it anti-symmetric. We therefore 
construct the spin-connection matrix as shown below: 


/ 


co = 


0 


2 m 


(r+m) 2 

) 


1- 


r+m 


<J- 


2 m 2 1 
(r+m) 2 

0 


m -a 3 


-(i -^y -(i- 


r+m 


+ 3 \ 


m a 3 


r+m 




a — 


r+m 


a 


r+m 

o 

O _ 2m 2 

\ (r+m) 2 




+^<T 2 

r+m 

_ O _ 2m 2 

V f ^ 


(r+m)^ 

0 


K 


(4.2.3) 


/ 


If we view the spin connections as a linear combination of self dual and anti-self dual tensors, 
then we can accordingly seperate out the self and anti-self dual components as con = ^. 

J L J L J 

To this end, we can split the spin connection matrix (4.2.3) into two separate components: the 
self dual and the anti-self dual parts 

; (+) = -^(o+i + +?72 + 0 + 3 ) = 

T 2 m 2 \ , ( 1 


co' 


CO( ) = 


(r + m) 2 


0+1 


m 


r + m 


(0+2 - +773) 


(4.2.4) 


For reference, we have the t’Hooft symbol matrices 7/+ which exhibit the su( 2) Lie algebra: 


[Vi,Vj] = ~^ij k dk 

The curvature tensor can be decomposed into self and anti-self dual parts itU = 4?+++■ \ 
where according to Cartan’s 2nd equation, R = dco + co A co. Thus, we can write the spin 
connections as as a linear combination of self and anti-self dual t’Hooft symbols giving us the 
self-dual and anti-self dual spin connections described in (4.2.4). Consequently, according to 
(4.3), the curvature tensor vanishes for the self-dual part: 

++) = d++) +++) A+ + ) = — ^ (dcd ++•£ cr 7 A cx fc ^77j = 0 (4.2.5) 

Only the anti-self dual curvature remains, reflecting the Taub-NUT’s anti-self dual nature. 
We make our job easier by writing the spin connection as = co[ ^ ^ 

w (_) = \ ( al Vi - a ' 2f l2 + v 3 m) + { - T-z ~+ ^ ~ + 773 ) ) (4.2.6) 

2 K [ (r + m) z r + m J 

where one can verify that co[ ^ will follow the same rule as in (4.2.5). 
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This allows us to compute the anti-self-dual curvature is given by 
2m 


= 


(r + m) 3 


rj i (e° A e 1 - e 2 A e 3 ) 


+ { - ^(e° A e 2 - e 3 A e 1 ) + %(e° A 

(r + m) 6 l ' v 


e 3 — e 1 A e 2 


(4.2.7) 


where we can see from the signs attached to the dual components that the curvature derived 
from Taub-NUT metric is clearly anti-self dual, as shown in [16]. This also lets us conclude 
that it is an instanton. To elaborate further, we can show that only S?7(2)_ gauge fields are 
embedded within the spin-connection components as shown below: 




"jlV 


A (±)» — \ 71 {F)i UJ 

A — ^'[iV ^ [IV 


(4.2.8) 


#)' = -- 


^(+) 2 = 


A(-) 1 = 


Al(-) 2 = - 


1 - 


4m 2 


r — mo 
r + m 2 


(r + m ) 2 
2 


o 


(4.2.9) 


^(+)3 = _ 


O' 


AH 3 = 


r — m o^ 


r + m 2 


while the field strengths are given by: 

r(-) = ri^ k F ^ 

'’[iv '/[iv k 

F 11 = Rqi = —i?23 = 

F ( )2 = i?02 = —^31 = — 


2m 


F 1 ' ' >:i = R 03 = —Rl2 = 


(r + m) 3 
m 

(r + m) 3 

m 

(r + m) 3 


7?(±)j _ _ TJ (±)i T> 

1 — a' i1u '’ ,w 

(e° A e 1 - e 2 A e 3 ) 


( e o Ae 2_ e 3 Ae i) 

( e 0 Ae 3_ e l Ae 2) 


(4.2.10) 


(4.2.11) 


where it is obvious that due to the absence of self-dual curvature, there are no SU{ 2) + 
gauge fields, ie. F = 0 and thus field strengths are anti-self dual (F = — * F) which 
of course, coincide with the curvature tensor (4.2.7). In terms of 2-forms, the independent 
components are given by : 


R (-) _ p(-) __p(-) _ 

-tVnni — i loQOQ — to 123 — 


2m 


0101 - “2323 - -■'0123 “ ^ + m y 


p( — ) m — ; pv — ; _ 

il non9 — 91 9 — -*Vl919 — 


>(-) 


(-) 


m 


"'0202 


9313 


0213 


(r + m) 3 


k>( _ ) — 7?( _ ) — — — 

Jinqnq — -tr-io-io — JUi-riQ — 


*'0303 


9212 


(-) _ 


m 


0213 


(4.2.12) 

(4.2.13) 

(4.2.14) 


(r + m) 3 

This lets us compute the Ricci tensors and scalar in accordance with the formula: 

R :i- g J, n, ikl S"-R ijkl R = S ik n ik (4.2.15) 


^00 — R-ll — R-22 — 77-33 — 0 R — 0 (4.2.16) 

Since Ricci tensors vanish, Taub-NUT is clearly a vaccum solution of Einstein’s equations. 
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4.3 Topological Invariants 

Topological invariants are analogous to an overall charge distributed in the manifold. In the 
gravity side, there are two topological invariants associated with the Atiyah-Patodi-Singer in¬ 
dex theorem for a four dimensional elliptic complex [17, 18]: the Euler characteristic x{M) and 
the Hirzebruch signature r(M), which can be expressed as integrals of four-manifold curvature. 


Recall that in electromagnetic theory, the Held action is given by: 

5 = ~~ [ dn Fijpi = f F A F 


where 


1 


F = -Fijdx 1 A dx J 


and 


e ljkl dQ, = dx l A dx j A dx k A dx l 


The equations of motion are obtained by solving for minimum variation of the electro¬ 
magnetic field action. We merely apply these equations to compute topological invariants as 
integrals analogous to action. We can write for the general lagrangian: 


£ = c abcd R ab A R cd = c abcd F£ )rn F£ )n rif )m V { ± )n £ ijkl dn 
= ±2 dnc abcd F^ )m d c A { ± )rn (where e ijkl dn = e i A e j A e k A e l ) (4.3.1) 

Applying Lagrange’s equation gives the contracted Bianchi identity for curvature as: 

(4 ' 3 ' 2) 

Conversely, we can say that Bianchi identity for SU (2)± gauge fields is at the root of the 
invariance of topological quantities. One can verify this starting from (4.3.2) and then working 
in reverse order to obtain the invariants. 

Given that the boundary integral vanishes, the overall invariant is computed only from the 
bulk part. For non-compact manifolds like Taub-NUT, there are additional boundary terms 
neither separated into self-dual nor anti-self-dual parts unlike the volume terms. They are the 
so-called eta-invariant rjsipM ), given for k self-dual gravitational instantons by [19] 


ris(dM) 


2e | (k — l)(k — 2) 
3 k + 3k 


e = 0; ALE boundary conditions 
e = 1; ALF boundary conditions 


(4.3.3) 


Since Taub-NUT is a ALF hyper-kahler four-manifold it has a non-vanishing eta-invariant 
which is equal to — |. According to calculations described in [20], upon applying curvature 
components of (4.2.7), the Euler characteristic x and the Hirzebruch signature complex r are: 

X( M ) = ^2 J e abcd R ab A Red = 1 (4.3.4) 


Tbulk (A7) — 


1 


12vr 2 


dd-ab A dd a b 


IM 


/ a<b 


2 

3 


(4.3.5) 


t(M) = T bu ik(M) + rjs(dM) = 0 

One could say that the general form of various topological invariants can be written as: 


c{m)= ^L 


c abcd R ab A R cd = 


M 


~j~2 Im F “b ( * F<lb ) ^ ^ ( Euler Char •) 

^ f F uF ab - 
JM rabr > 


C abcd _ gacgbd (Hirzebruch Sign.) 


where c abcd is contracting tensor defined in respect to the relevant circumstances. 


(4.3.6) 
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5 Killing-Yano tensors and the Taub-NUT metric 

There are tensors quadratic in momenta and conserved along geodesics, expressed as a vector 
K whose components transform among themselves under 3-dimensional rotations. They are 
very similar to the Runge-Lenz vector in the Kepler problem with components: 

K® = \K^p^ v (5.1) 

Provided that J° / 0, such vectors usually satisfy the following property: 



where if (J°, J, H, K) are all constant, the 3-dimensional position vector r lies in a plane. 

V. J 

Using (5.2) and the relation <7° =-, we can see that: 

r.K = T rH + ^J 2 - (J 0 ) 2 ) (5.3) 

In Taub-NUT geometry, there are also 4 completely antisymmetric Killing tensors known 
as Killing-Yano (KY) tensors. Three of these are complex structures, realizing quaternionic 
algebra since the Taub-NUT manifold is hyper-Kahler. The fourth is a scalar with a non¬ 
vanishing held strength and it exists by virtue of the metric being of Petrov type D. Their 
existence is implied by a triplet of symmetric 2nd rank Killing tensors called the Stackel-Killing 
tensor satisfying: 

0(A K®, = 0 (5.4) 

We will examine properties of KY tensors relevant for studying Taub-NUT symmetries. 
Before that, let us list some references that initiated the study of such dynamical symmetries. 

Dynamical symmetries of the Kaluza Klein monopole were discussed in detail by Feher 
in [21]. The dynamics of two non-relativistic BPS monopoles was described using Atiyah- 
Hitchin metric (Taub-NUT being a special case), the corresponding 0(4)/0(3,l) symmetry 
discovered in [22], and applied to calculate the underlined motion group-theoretically in [23]. 
The symmetry was then extended to 0(4,2) in [24] and [25]. In [24] Gibbons et. al discussed 
dynamical symmetries of multi-centre metrics and applied the results to the scattering of BPS 
monopoles and fluctuations around them, giving a detailed account of the hidden symmetries of 
the Taub-NUT. The hidden symmetries in large-distance interactions between BPS monopoles 
and of the fluctuations around them are traced to the existence of a KY tensor on the self-dual 
Taub-NUT. The global action on classical phase space of these symmetries was discussed in 
[26] and the quantum picture involving the ’’dynamical groups” 50(4), 50(4,1) and 50(4, 2) 
was also given. A comprehensive review of the dynamical symmetry can be found in [27]. 
Supersymmetry and extension to spin has also been studied in [28, 29]. 

5.1 Yano and Stackel tensors 

We can construct these KY tensors in terms of simpler objects known as Yano tensors that 
are antisymmetric rank 2 tensors satisfying the Killing like equation. Thus, the covariant 
derivative is antisymmetric over permutations of all possible pairs of indices. This allows us 
to write the covariant derivative of the Yano tensor in terms of the cyclic permutations as: 

Uu = V + V = o (5.1.1) 

V m /„a = Vvfx/j, = VxUv = V [ij,fvx] = A + V^/a m + Va fnv) (5.1.2) 
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We can also construct symmetric Killing tensors of rank 2 by symmetrized multiplication: 



These symmetric Killing tensors satisfy the KY condition (5.4). The Taub-NUT manifold 
admits 4 such KY tensors, given by a scalar /° and three components that transform as a 
vector p V i = 1,2,3. We can form triplets of symmetric Killing tensors as in (5.1.3), given 
by setting a = 0 and b = i: 

K® = Kff = I (/° A /L + f^fl) * = 1,2,3 (5.1.4) 

Using (5.1.1) we can see how they obey (5.4) as follows: 

= 0 

V(rC)=0 ==• V b K U = = 0 (5-1-5) 

Thus, we can feel assured that (5.4) is satisfied by this symmetric Killing tensor. This 
allows us to construct the tensors of (5.4) that are quadratic in momenta, showing how to 
get Stackel tensors from the KY tensors. However, since the KY tensor is anti-symmetric, it 
cannot be used to form polynomials with components of the same vector. Thus, it will have to 
be a mixed product of components of different vector quantities, as found in case of the angular 
momentum, a product between one position and one momentum component each. Applying 
Holten’s algorithm yields the Killing equation in (5.1.1). 


5.2 Euclidean Taub-NUT 

The Taub-NUT metric [30] admits four such Yano tensors written as the following 2-forms: 

f° = 4 {dip + cos 6 d<p) A dr + 2 r(r ± l) (r ± 2) sin 0 d6 A dp (5.2.1) 

P = ±4(ch/’ + cos 9 dp) A dx l — e l j k f{r) dx^ A dx k , Vi, j, k = 1, 2,3 (5.2.2) 

One can always find Killing tensors embedded within conserved quantities, as evident from 
the Poisson Brackets of any conserved quantity expanded ala Holten algorithm. The coefficient 
from Laplace-Runge-Lenz vector is analogous to the Killing-Stackel tensor Kij, so we can argue: 

Q< 2 > = KijlTIP = ^cj;fu l U j (5.2.3) 

Now the angular momentum co-efficients according to (2.2.10) are: 


C (0) = q g jk (x) -9 k = -g im (x)e m jk 6 k xi 


(5.2.4) 


If we write cP' 1 = fi k 0 k (see section 8.1 in APPENDIX), using Holten’s Algorithm gives: 

Vj-cf } = vp lk e k = -g im {x)£ m jk e k 
VU, (1) + VjCV = 0 =A (vj jk + Vj/ifc) e k = 0 


which is the Killing equation (5.1.1). Thus, we can say that the KY tensor is 

fjk = 9jk(x) => /ofe = 

fjk = £l jk => P = £ l j k e J A e k 
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(5.2.5) 

(5.2.6) 


such that the square of it gives the Stackel tensor 


(5.2.7) 

This shows how Killing tensors are embedded within the conserved quantities. We can 
choose four combinations of three indices out of the available four. Since Taub-NUT can be 
written in an alternate form given by (2.1.1), the vierbeins of the metric are given by: 

= e , = ^ ixi ( 5 . 2 . 8 ) 

vW) 

So, according to our theory, we should have 


f = -e jk e i A e k + e° A e k 

= — e L jkf{r)dx^ A dx k ± 4 [dip + A.dxj A dx 1 


(5.2.9) 


This result so far is comparable with the result (5.2.2), so we have a possible method for 
constructing Killing-Yano tensors from the coefficients of conserved quantities. Their covariant 
exterior derivatives and their properties are given by: 

Df° = V 7 /°„dz 7 A dx 11 A dx v = r(r ± 2) sin 6dr A d6 A dcj) (5.2.10) 

Df = 0, Vi = 1,2,3 (5.2.11) 


From the results above, we can infer that the covariant derivatives hold following properties: 


V 7 /J„ = V„/° 7 = V 7 4„ = 0 * = 1,2,3 (5.2.12) 

showing that they obey the condition for covariant derivatives of KY tensors. As shown 
in (5.1.4), these tensors can form a symmetric triplet or a vector of Killing tensors. They also 
exhibit the mutual anti-commutation property: 


= + = -2<W 

[fj j ] =f i f j ~f j f i = 2e i hf k 


(5.2.13) 


proving that they are complex structures realizing the quaternion algebra. This implies that 
the 2-forms f l are objects in the quaternionic geometry and possibly hyperkahler structures. 
This leads us to the next section where we examine the hyperkahler structure of the Taub-NUT. 


5.3 Graded Lie-algebra via Schouten-Nijenhuis Brackets 

We will now see if the KY tensors of the Taub-NUT metric exhibit Lie algebra under the 
action of Schouten-Nijenhuis Brackets. If they do, it would allow us to form higher order KY 
tensors from lower order ones of rank greater than 1. In particular it is noteworthy in this 
context that, Kastor et. al already found that KY tensors on constant curvature spacetimes 
do form Lie algebras with respect to the SN bracket [31]. 


The Schouten-Nijenhuis Bracket (SNB) is a bracket operation between multivector fields. 
The SNB for two such fields A = A tli2 - tm AfcLi ; B = d Jk . is given by 


£iai...a m+7l _i _ ^ ai-..«m+n-l 

= mA C ^ ai '" a " l ~ 1 'V ,B am "' arn+n ~ -|- 7l( — -Q mn ^c[ai...a n -i^7 J^a n ...a,m+n- 1] 


(5.3.1) 


This new tensor is completely antisymmetric, fulfilling the first requirement to be considered 
a KY tensor. All that remains is for its covariant derivative to exhibit the same Killing 
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equation (5.1.2) relevant to such tensors. Now, we will use an important identity (see (8.3.1) 
in APPENDIX) for KY tensors: 

VaV^,..^ = (-l)^^R [bci{ /K C2C3 ,,, Cn]d (5.3.2) 

we get upon applying to the covariant derivative of this new tensor: 

VbCa 1 ...a m+n _ 1 = — (m + n) (V c A[fe ai ...a m _ 1 )V B am a rn+n _ 1 ] 

c d (5.3.3) 

— (m + n)A [ ai ...a m _i-^|6d|ca m -®a m+ i...o m+n _i] 

The 1st term easily shows anti-symmetry of index b with other indices, but the 2nd term 
exhibits it only under certain circumstances. One could say that by symmetry properties of 
the curvature tensor, in maximally symmetric spaces it could be expressed as: 

Rabcd(x) = f(x)gij(x)e\ b £ j cd = f(x){g ac (x)g bd (x) - g a d(x)g bc (x)} (5.3.4) 

So, for cases of constant curvature f(x) = k , we could write 


{Rabcd) const = k{g ac {x)g bd {x) - g a d{x)g bc {x)} (5.3.5) 

Thus, upon applying the constant curvature formula of (5.3.5) to (5.3.3), we will get 
VbCa^a^ — — (m + n) [(V c A^ ai Clm _ 1 ) V -B am ...a m+n _i] 


- kAi 


[<Jl R CLjn...dm-hn — lb 


= V \ b C, 


(5.3.6) 


[6 CL±...CLm-\-n — l\ 


Clearly this matches the property eq.(5.1.2) expresses, showing that it is also a KY tensor. 
So the SNB of any two KY tensors in spaces of constant curvature is also a KY tensor. 

However, as evident from (4.2.7), the curvature of the Taub-NUT metric is not constant, 
allowing us to conclude that its KY tensors do not exhibit Lie algebra under SN Brackets. 
Thus, we cannot produce higher order KY tensors using the lower order ones for the Taub- 
NUT as shown in [32], So, we are limited to the set of four available rank two KY tensors. 


6 Hyperkahler structure and the KY tensors 

Now we will consider the hyperkahler structures related to the Taub-NUT metric. To begin 
with, we will define both, kahler and hyperkahler structures. 

Kahler manifold: If a complex manifold M has a hermitian metric g and a fundamental 
2-form oj which is closed (dco = 0) then M is a Kahler manifold and w is a Kahler form. 

The connection between the metric g and the Kahler form u is: 

= (Jg)^ ( 6 - 1 ) 

where J is the complex structure, for which J 2 = —1. 

Hyperkahler manifold: If M is a hyper-complex manifold with a hyper-Hermitian metric 
g and a triplet of fundamental forms oj which are closed (did = 0) then M is a Hyperkahler 
manifold. It is the same as the Kahler manifold except that there are more than one type of 
complex structures. In case of 4 dimensions, there are 3 such integrable complex structures 
( i,j,k ), and they obey the algebraic relations: 

i 2 = j 2 = k 2 = ijk = —1 (6.2) 
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This would also imply that there are corresponding number of different 2-forms available 
in this case, known as the Hypei'kahler forms, given by: 

< = = (6-3) 

where g\ v is the hyper-hermitian metric and J l ^ is the almost complex structure exhibiting 
quaternion algebra 

Jry ,Jft 8(yfjl T £afP J 7 (6.4) 

and thus, we can see that the hyperkahler structures exhibit the same algebra: 

[A ]% v = 2^*4- 

(«w) = <W. = (4"-sn)s'’ i (4 °-= jV4°h™ = 

= ([A •»*]») = 2(«V*»)„ = 2^^ 


These complex structures originate from the t’Hooft symbols which have 3 self dual and 3 
anti-self dual components. That means we could have six different symplectic 2-forms. The 
almost complex structures J l can be represented by t’Hooft symbols, which themselves can be 
given by linear combinations of antisymmetric tensor e l j k and delta function 5 l j. 


n _ j 

J jk — 8 jk 


±\{5° j 8 i k -5° k 8 i j ) 


Thus, we can argue that hyper-kahler structures given by (6.3) are: 

Ujk = ( J l g) jk = 9j n(x) e in k ± l -(8* n 8\ - 8° k 8 in ) 

As introduced in (2.1.2) and following [33] we shall take a different form of the Taub-NUT 
ds 2 = V(r) Sij dx l dxi + U _1 (r)(dr + a.dr) 2 
for which, the vierbeins, in a similar fashion to (5.2.8) are given by 


(6.5) 

( 6 . 6 ) 

■NUT 

(6.7) 


e° = 


4 (dr + a.dr) 


e* = y/V(r) dx 1 


Thus, remembering that g = 8^ e* (g) e 3 the hyper-kahler forms are given by: 

u/ = w* fc dx J A dx k = J l j k e •? A e k 


( 6 . 8 ) 


uj 1 = 


£ i jk ±-(8 0 j 8 i k -8° k 8 i j ) e j Ae k = e i jk V(r)dx j A dx k - e° A e* 

uj 1 = e l j k V (r)dx^ A dx k ± (dr A dx l + a n .dx n A dx l ) (6-9) 

For the Taub-NUT, choosing only anti-self-dual components for V(r) = l + — and restricting 

r 

<7 to lie on a plane (<r = ( 0 , <72, <73)), the reduced symplectic forms are: 

c j 1 = dx l A dr + c^dx 1 A dx 2 + a^dx 1 A dx 3 + H—^ dx 2 A dx 3 

a; 2 = dx 2 A dr + a^dx 2 A dx 3 — + -^) dx 1 A dx 3 


( 6 . 10 ) 


a; 3 = dx 3 A dr — c^dx 2 A dx 3 + ^ + -^ dx 1 A dx 2 
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This construction of hyperkahler structures is similar to how spatial KY tensors were 
deduced, proving that the KY tensors are the hyperkahler structures of the Taub-NUT metric. 


Comparison between Killing-Yano Tensors and Hyperkahler Structures 


i 

Killing-Yano tensor /* 

Hyperkahler structure uA 

i 

±4 (dil: + A n dx n ) A dx l — e l jkf(r)dx J A dx k 

±(dr + a n .dx n ) A dx 1 + £ l jkV(r)dx^ A dx k 

1 

T^dx 1 A (dip + A n dx n ) + ^1 + -^ dx 2 A dx 3 

dx 1 A (dr + a 2 dx 2 + cr^dx 3 ) + dx 2 A dx 3 

2 

q=4 dx 2 A (dip + A n dx n ) — ^1 + dx 1 A dx 3 

dx 2 A (dr + a^dx 3 ) — 4—^ dx 1 A dx 3 

3 

q=4 dx 3 A (dip + A n dx n ) + ^1 + dx 1 A dx 2 

dx 3 A (dr — a 2 dx 3 ) + (l -\—^ dx 1 A dx 2 


Few points are worth mentioning here. By studying the G 2 holonomy equation for biaxial 
anti-self dual Bianchi IX base Gibbons et.al [34] found that the associated first order equa¬ 
tions satisfied by the metric coefficients yield the self-dual Ricci flat Taub-NUT metrics where 
50(3) C U( 2) rotates the three hyperkahler forms as a triplet. 

7 Discussion 

In this article we see that the Taub-NUT is comparable to Euclideanized Bertrand spacetime 
with magnetic fields due to the shared geometry and conserved quantities, and a dual config¬ 
uration as either Oscillator or Kepler systems. Identical conserved quantities are connected 
to identical symmetries and Killing tensors embedded within. These are the Killing-Stackel 
and Killing-Yano tensors embedded as co-efficients within the Laplace-Runge-Lenz and an¬ 
gular momentum vectors respectively. The Killing-Yano tensors exhibit quaternionic algebra, 
hinting at a link between them and hyperkahler structures matching the form of the KY ten¬ 
sors derived from the angular momentum. This confirms that the KY tensor and hyperkahler 
structures are the same for Taub-NUT. Since symmetries of a spacetime are unaffected upon 
euclideanization, we can expect that all properties arising from shared symmetries are also 
exhibited by Bertrand spacetimes with magnetic fields. 

Taub-NUT is a special case of the anti-self-dual Bianchi-IX spaces [14], derived by solving 
the dynamical equations that emerge from applying the settings for this case to the related 
classical Darboux-Halphen system. The shared geometric properties, including the Ricci flow, 
integrability aspects and integrable reductions to Painleve systems can be explored to some ex¬ 
tent. In special situations, self-dual Einstein Bianchi-IX metrics reduce to Taub NUT de Sitter 
metric with two parameters of the biaxial solutions respectively identified as the NUT param¬ 
eter and the cosmological constant. Computing its curvature confirms that the Taub-NUT 
is anti-self dual, reflecting its instantonic nature, and expectedly, Ricci-flat with topological 
invariants to compare with other possible diffeomorphically equivalent Ricci-flat manifolds. 
According to Kronheimer classifications [35, 36] all hyperkahler metrics like Taub-NUT in 
four dimensions are always anti-self dual, so the hyperkahler quotient construction, due to 
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Hitchin, Karlhede, Lindstrom and Rocek [37] carries an anti-self dual conformal structure, 
allowing Penrose’s Twistor theory [38] techniques to be applied in this case. 

Recently works in emergent gravity [39] aim at constructing a Riemannian geometry from 
17(1) gauge fields on a noncommutative spacetime. This construction is invertible to find 
corresponding 17(1) gauge fields on a (generalized) Poisson manifold given a metric ( M,g ). 
There are already detailed tests [40] of the emergent gravity picture with explicit solutions in 
both gravity and gauge theory sides. Symplectic U(l) gauge fields have been derived starting 
from the Eguchi-Hanson metric in four-dimensional Euclidean gravity. The result precisely 
reproduces the U(l) gauge fields of the Nekrasov-Schwarz instanton previously derived from 
the top-down approach. To clarify the role of noncommutative spacetime, the prescription has 
been inverted and Braden-Nekrasov 17(1) instanton defined in ordinary commutative space- 
time was used to derive a corresponding gravitational metric just to show that the Kahler 
manifold determined by the Braden-Nekrasov instanton exhibits a spacetime singularity while 
the Nekrasov-Schwarz instanton gives rise to a regular geometry-in the form of Eguchi-Hanson 
space. This result implies the important role noncommutativity of spacetime plays in resolv¬ 
ing spacetime singularities [41] in general relativity. Some relevant studies related to emergent 
nature of Schwarzschild spacetime was also performed in [42], 

One may wonder if we can get 17(1) gauge fields in the same way from the Taub-NUT metric. 
A critical difference from the Eguchi-Hanson metric [17] is that the Taub-NUT metric (2.1.1) 
is locally asymptotic at infinity to R 3 x S 1 , and so it belongs to the class of asymptotically 
locally flat (ALF) spaces. Thus, the Hopf coordinates cannot represent the Taub-NUT metric, 
and it is difficult to naively generalize the same construction to ALF spaces. From gauge 
theory perpective, it may be related to the fact that ALF spaces arise from NC monopoles [43] 
whose underlying equation is defined by an S^compactification of the self(anti)-dual-instanton 
equation, the so-called Nahrn equation. We will discuss in [44] a possible generalization to 
include the Taub-NUT in the bottom-up approach of emergent gravity. 

It is known that only for special choice of the NUT parameter we get a regular metric, but 
generally, one encounters singularities at either end point of the 4-dim radial coordinate. In the 
most generic case, although with a particular choice of the period for the azimuthal angle, one 
can get away with the bolt-singularity, the NUT singularity (co-dimension 4 orbifold singular¬ 
ity) stays, possibly admitting an M theory interpretation associated with the corresponding 
non-abelian gauge symmetries [45]. 

Recently, Ricci flat metrics of ultrahyperbolic signature were constructed [46] with 1- 
conformal Galilei symmetry, involving an AdS 2 part reminiscent of the near horizon geoemtry 
of extremal black holes. Similarly, it should be interesting to see if Taub-NUT spaces are as- 
sociable with geodesics that can describe second order dynamical systems. Perhaps the most 
interesting issue will be to explore whether something like “Taub-NUT/CFT” correspondence 
can be conjectured. 
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8 APPENDIX 

Important computations and derivations of this article are provided in this section. 


8.1 Basic Killing tensors from Holten’s Algorithm 

Angular Momentum 

(n) 

If we choose to set Cl. 7 =0, V n > 2, we get the Killing equations: 

V^cjf =0 (8.1.1) 

There are two parts of this solution we shall study in detail. We can write (8.1.1) as: 

ViCf } + VjCf : ’ = 0 => ViCf ] = -V j c[ l) (8.1.2) 

This is an anti-symmetric matrix, written as = —6ji- Further elaboration gives: 

9ij(x) = £ijk{x)0 k = gim(x)E m jk e k 
-VjC^ =g im (x)£ m jk e k => cf } = -g im (x)e m jk 0 k X j 
Thus, we have the rotation operator as the 1st order co-efficient: 

Cf } = -g im (x)e m jk 6 k xi (8.1.3) 

Applying this co-efficient into the 1st term of the power series, we get: 

Q {1) = cf } n l = - gim {x)£ m jk e k x^ 

=» L.9 = ~(e ijk U i x j )e k = (xxU).e 

L = x x II (8.1.4) 

This eventually becomes the conserved quantity known as the angular momentum. 


Laplace-Runge-Lenz vector 

(n) 

Now when we choose to set Cj,.| =0, V n > 3, we get the Killing equations: 

7 iOS + V t c<, 2) = 0 


ViC]*' + V,ci 


- ( 8 - L5 ) 
as we can see, (8.1.5) perfectly matches the property of the Killing Yano and Killing Stackel 
tensors. The Runge-Lenz like quantity is given by a symmetric sum as shown below: 


A x [B x C) = £u rn £ rn j k A l B j C k £ii m £ rn jk = 8ij5ik - 5 ik 5ij 


( 8 . 1 . 6 ) 


VfcCf } = £u m {x)£ m jk{x)n l + (ie j) 

= (2 gij(x)g k i(x) - g ik (x)gji(x) - gu(x)g kj (x))n l x k 
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( 8 . 1 . 7 ) 


C\f = ( 2gij(x)n k - gik(x)rij - g k j{x)rii)x k 


As before, applying this co-efficient to the 2nd order term in the power series gives 

12 . 


^ (2) h - 

= N.n = {|n| 2 ® - (n.®)n}.n = {n x (x x n)}.n 


Q (2) = -qf WW = {|n| 2 (n.a:) - (U.x)(U.n)} 


n = n x {x x n) 


( 8 . 1 . 8 ) 


This quantity is a term that is present in another conserved quantity known as the Laplace- 
Runge-Lenz vector. Having found the two familiar types of conserved quantities, we can now 
proceed to see what it looks like for the Taub-NUT metric. 


8.2 The Bohlin transformation 

The Bohlin transformation that maps the co-ordinate system on a plane, is given by: 

/:*—>£“ = (z a ) 2 = Re i * => z = £* 


( 8 . 2 . 1 ) 


Now we must note that another Noether invariant, the angular momentum will change form 
under this transformation. We re-parametrize to preserve the form of angular momentum. 


i = r z e = \z\ 2 e = if | V 


|2 -/ 






The velocity and acceleration can be given as: 


z a = v^ T r , = ^(r) l r' 


z = 


2 (£a): 

l|£i|{(r)V} = l^rr" + l(r) J i«'i 


( 8 . 2 . 2 ) 

(8.2.3) 

(8.2.4) 

(8.2.5) 


The equation of motion for a Harmonic Oscillator eventually becomes: 

4^ 


m < i-l^ r f a,, + -ACY 2 \Z '\ 2 \ = -k{e) 


(£“)• 


=* \z\ 2 *; a,, + h a \?\ 2 = --t a => a a " = -(l\e\ 2 + 

2 m \2 m J \t,\ z 


( 8 . 2 . 6 ) 


The Hamiltonian TL of the oscillator can be re-written to complete the transformation: 


m,., 0 k. l9 mf 1 ..,., 2 k\... 

« = yN 2 +2 |2|2 = T U | «'l 2 + -)!£! 


, |£"'| 2 , 2*\ = 





+ 


m J | 2 


l?l 5 


|f'| 2 2k\ _ m J_ _ J_ 

2 m) m |f| K |f| 

(8.2.7) 


showing that it restores the central force nature of the system, giving us the equation of 
motion for inverse square law forces. 
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8.3 Double derivative of Killing-Yano tensors 

Similar to Killing vectors, rank n Killing-Yano tensors exhibit a curvature equation 


(V a V 6 - V(,V a )-K Cl ...c„ = '^2 Rabci d Kci...d...c 


i =1 


For the LHS, by permuting the indices according to the rules, we will get 
LHS = (V a V fe V 6 V a )K Cl ... Cn = — V a V Cl -fQ, C2 ... Cn + VbV Cl -K aC2 ...c n 

n 

2V C1 Vf,K aC2 ...c n Rabci K-dc2...Cn ^ ^ (-^6ciCi K a c2...d...Cn Rac±Ci K-bc2...d...Cn ) 


i=2 


Rabci ^dc2...c n Raba d K C1 ...d... C 


i=2 


2V Cl VfcA aC2 ... C71 — 2 R a bc\ Kdc 2 ...c n + ^ ^ ( Raaci K b c 2 ...d... Cn Rbcia R-ac 2 -..d...Cn "P Rabci A-ci...d...c n ) 

z —* V--------■ V---' 


i=2 


II 


V cl V 6 JF aC2 ...c„ = V cl V[ fe A: aC2 ... C7i] , W C1 := V C1 V [bK aC2 ... c „]e“ A e 6 A e C2 ... A e c 
With I and II above, we can create the 3-forms 

I R abCi d e a A e 6 A e c ‘ = - ( Robe* + RbcJ + R Ci ab d )e a Ae b Ae Cl = 0 

II ->■ Rac 1 c d K b c 2 ...d...c ri e a A e b A e Ci = -R Cl ac d Kbc 2 ...d...c n e- a Ae b A e Ci 

® A A zd*-'* 


= —Rciab K a c 2 ... d ... Cri e A e A e 1 


Thus, II will become: 


r, {RciCib d Kac 2 ...d...c n + Rbc 1 c d Kac 2 ...d...c n )e a A e b A e Ci = ^ Rc i bc- l d Kac 2 ...d...c n e a a e b A e Ci 

i=2 
n 

= J2 R abc 1 dK dc 2 ...c i ...c n e a A e fe A e Ci = (n - 1 )R a bc 1 d K dc2 ... Ci ... Cn e a A e b A e Ci 


i=2 


i=2 


Applying this result back in the main equation, we get: 

2V C1 V b K aC2 ...c n e a Ae b A e Ci = (2 R abci d K dc2 ... Cn + (n - l)iU cl d -FQ C2 ... c Je a A e 6 A e Ci 
=A 2V C1 V b K aC2 ...c n e a A e b A e Ci = (n + l)R a bc 1 d Kdc 2 ...c ri e a A e b A e Ci 

Finally, we get the double-derivative of KY tensors as: 


N a N h K clC2 ... Cn = (-l)-+i^±li? [ftci|a| ^ C2C3 Cn]d 


(8.3.1) 
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